Creation of maximally entangled photon-number states using optical fiber multiports 
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We theoretically demonstrate a method for producing the maximally path-entangled state 
(1/^/2) (|A*', 0) + exp[iA'"(;!>]|0, A'^)) using intensity-symmetric multiport beamsplitters, single photon 
inputs, and either photon-counting postselection or conditional measurement. The use of postse- 
lection enables successful implementation with non-unit efficiency detectors. We also demonstrate 
how to make the same state more conveniently by replacing one of the single photon inputs by a 
coherent state. 

PACS numbers; 03.67.Mn, 42.50.Dv, 42.50.St, 06.30.Bp 



There is great interest in the development of a range of 
nonclassical optical states for quantum technology appli- 
cations jlj , including applications in quantum computing 
0; quantum cryptography [3; quantum metrology and 
lithography |3. Isl IRT^ IslM llo| : and quantum imaging 
[lll| . A development of particular interest is the "quan- 
tum Rosetta stone" , which elucidates the isomorphism 
between Mach-Zehnder interferometers (MZIs), Ramsey 
spectroscopy and quantum phase gates. It has been 
shown theoretically that a MZI produces interferometric 
fringes with a period A/TV when employed with a maxi- 
mally path-entangled iV-photon state, 
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within the interferometer, where the ordered pair in the 
kets represents the photon number in each of the two 
interferometer arms. The period is therefore 1/N times 
smaller than that for a MZI employed with a coherent 
state or single photon input of wavelength A 0- This 
has been experimentally verified for iV = 2 and an 
equivalent result obtained in Ramsey spectroscopy [l^ . 
These reduced-period fringes can be used to measure a 
phase shift in one arm of the interferometer with a phase 
uncertainty A(j> — 1/N at the Heisenberg limit for all 
N and (j) 0- The same quantum state |'0)noon 
been suggested for use in generating spatial interfero- 
metric fringes with a period TV times smaller than single 
photon fringes, with potential application to beating the 
diffraction limit in lithography [3. The reduced period 
of a spatial interference pattern has been demonstrated 
for the case N = 2 There is significant interest in 

the extension of these techniques to larger iV, and there 
exist a number of proposals for achieving t his, each of 
which is experimentally intensive [TsL [T3. Iisl [la | . 

In this paper, we theoretically demonstrate a method 
of producing the quantum state |V')noon' ^'^^ general iV, 
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FIG. 1: (a) Conceptual diagram of 2A'^-port intensity- 
symmetric beamsplitter, with input and output modes la- 
beled, (b) Mach-Zehnder interferometer, for measuring phase 
shift (f), created with a 2A^-port splitter and a 50/50 beam- 
splitter. Postselection occurs when the two detectors count 
a total of A'' photons between output modes and 1. The 
postselection guarantees that the state at B was the maxi- 
mally number entangled state, and fringes with a period 1/7V 
times the single photon wavelength are observed in N-photon 
coincidence. 



by means of a 2iV-port {N inputs, N outputs) beamsplit- 
ter [TtLITM Ii^ with symmetric splitting ratios, and pho- 
ton counting. The technique we present is less intensive, 
in that it uses a single beamsplitting device instead of a 
cascaded array of beamsplitter and measurement stages, 
and there is correspondingly a significant reduction in 
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the number of optical elements required. Further, it has 
the advantage of operating successfully with non-unit ef- 
ficiency detectors because it employs postselection. An 
intensity-symmetric 2iV-port beamsplitter is readily im- 
plemented in practice with a N x N fiber coupler, which 
are available commercially for various values of N up to 
TV = 32. A further development is the simplification of 
input state preparation by replacing one of the single 
photon inputs with a coherent state. 

The 2 A^-port itself is represented in Fig. QJa) , where 
it is employed with A^ single photon inputs. The input 
state can be written: 
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To obtain the output state, we act on the input state 
with the 2A^-port transformation operator, and use its 
unitarity and the time reversal symmetry of the network 
to write: 
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where a^(6j, . . . , bl^_i), abbreviated a[.(6^) henceforth, is 
the input operator for the fcth mode written in terms of 
the output operators In the case of the canonical 
2Af-port, the &^ are related to the input bosonic mode 
operators at by the matrix equation. 



at(6t) = M, 
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where the canonical 2A^-port matrix elements are given 
by the discrete Fourier transform jl8| |. 
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and fc, I vary between and A^ — 1. As an example, the 
canonical 6-port matrix is 



1 



1 1 



1 



V'i \ I g47ri/3 g27ri/3 

and the output state, with single photon inputs, is 
IV'3)b = + + 36^5165) |0,0,0) 



(6) 
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This state is itself highly entangled. It is also of a form 
that is somewhat suggestive of the IV') noon state. In 
general, to obtain the |V')noon state from that of Eq. JSJ, 
we invoke postselection using number-resolving photon 
counters ,2Qj|, as shown in Fig. ^b). The postselection 
occurs when the counters together count a total of A^ 
photons simultaneously. Due to the requirement to have 
a total of A^ photons in output modes and 1, the only 
terms retained in the expansion of Eq. Q are those with 
operator products of the form b\'^~'^b^i , with < fc < A^. 
It is important to note that the counting of A^ photons 
must occur after the final beamsplitter, otherwise the 
states |Af, 0) and |0, A") are distinguishable, because the 
"which-path" information is accessible. In principle, it is 
possible to replace the MZI with other experiments that 
utilize I "0) NOON I ^s long as the path information is not 
accessible. 

The postselection, performed as shown in Fig. Q^b), 
ensures that we know with certainty that whenever the 
detectors count a total of A^ photons in modes cq and ci , 
there are no photons in modes C2, . . . , cm-i and the state 
at B, the output of the 2Af-port itself, is known to have 
been 
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since the coefficients of the cross terms between &J and 
b\ vanish, as proved in the Appendix. The postselection 
ensures that this technique works, even with non-unit 
efficiency number resolving counters, since at most A^ 
counts are possible, and all counts less than A^ are re- 
jected. The state \'4')b^ is normalized by the coefficient 
K = N^/^/{y^{Nl)V2), and using Stirling's approxima- 
tion in the case of large A^ yields 1/K « V T:Ne~^ . The 
latter gives the scaling law for the production of the max- 
imally path-entangled A^-photon state using our method. 
This exponential scaling law is similar to the scaling laws 
for other proposed methods 0, Q, ^| of generating the 
state I V-) NOON- 

Our technique can also be employed to prepare 
|'(/')noon conditionally using projective measurement as 
opposed to postselection. The condition of A' photons 
in the output modes bQ,bi of the multiport is enforced 
by measuring vacuum on each of modes &2, ^7V-i, in- 
stead of counting photons at the output of the MZI. Con- 
ditioning on the measurement of zero photons has the 
obvious disadvantage of requiring unit efficiency detec- 
tors, but with the advantage of removing the destructive 
measurement on modes and 1. 

The non-postselected output state of the complete MZI 
can be calculated using the time-reversed creation oper- 
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ators as before, 
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but now the unitary operator is the transfer operator 
for the entire system, represented by the N x N matrix 
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$^M4g^, composed of the transfer matricies of 
the different elements of the optical network. The matrix 
$ represents the phase shift in Fig. ^b) and is the 
identity except for the first diagonal element which is e*"^ . 
The matrix M4 ox represents the final beamsplitter in the 
MZI and is the identity except for the four elements in 
the top left corner which are the matrix elements for the 
canonical 4-port, as per Eq. ©. Using the same proof 
as before, the (unnormalized) postselected state is 
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The form of this state agrees with the output state pre- 
dicted for a Ramsey ^spectrometer or MZI employed with 
the state |?A)noon -J^- Following Bollinger et al. 0, the 
phase measurement is made by measuring the parity of 
the count in the second detector, (— , which are the 
eigenvalues of the observable O: 

O(aS±al)^|0,0) = (aSTaI)^|0,0), (11) 

and which is equivalent to applying the Pauli operator 
to each of the N photons. The expectation value of O 
is then {6)^^ = (-1)^ cos(iV(?!)). Clearly the period of 
the interference is reduced by a factor of 1/ A'^ over single 
photon interferometric measurement of a phase shift 0, 
as expected for the maximally path-entangled state. It is 
straightforward to calculate the variance and deduce that 
the phase uncertainty is 1/N 0, which is at the Heisen- 
berg limit for the subensemble of postselected events. 

We now consider the application of the multiport fiber 
beamsplitter in the case where one input state is replaced 
by a coherent state. Using a modified version of Eq. (O 
with one Fock state replaced by a displacement operator, 
the output state of the 2A^-port at point B, in the absence 
of postselection, is 
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where Dk{a) is the coherent displacement operator act- 
ing on mode k. Postselection by counting N photons in 
the first two output modes, as before, produces the state 
IV')noon limit of small a, i.e. where is much 

smaller than all the other photon input rates in the sys- 
tem. To see this, observe that the Fock inputs can at 
most contribute A'^ — 1 photons to the first two outputs. 



In general, the — 1 single photon inputs contribute 
N - I- q (with Q < q < N - I) photons to the first 
two output modes, and q photons distributed somehow 
amongst the other output modes. The coherent state in- 
put will be required to contribute q + \ photons to the 
first two output modes to make a total of N photons in 
those modes, and trigger the postselection. The largest 
amplitude term(s) of the state with TV photons in the 
first two modes goes as a'^^^ . In the small a limit, we 
discard terms with q > Q and the product of coherent 
state operators in Eq. l(T^ becomes 
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The other terms of order a, i.e. a 6j (1 < p < iV — 1) are 
lost in the postselection because they do not contribute 
any photons to the first two output modes. With the 
approximation l|13|) . Eq. (|12|l collapses to ||SJ and the 
maximally entangled state is obtained. 

There exists an exact method (that does not constrain 
a) for realizing this same transformation |1, 1, . . . , a) — > 
I V') NOON' involving the use of photodetectors on all other 
outputs of the network of Fig. ^b), as well as the first 
two. When each of the detectors on output ports 2 to 
iV — 1 registers no count, simultaneously with the usual 
counting of N photons in the first two output modes, 
the approximate relation H13|l becomes exact. This is 
because all terms of order 2 or higher in a create at least 
one photon in an output mode other than or 1, and so 
are postselected away. 

An interesting research problem is the replacement of 
more than one of the Fock state inputs with coherent 
states. This is potentially a rich field of investigation, 
because of many potential phase and amplitude relation- 
ships that could be chosen between coherent input states. 
An open question is the effect of non-unit efficiency num- 
ber resolving detectors when employing our technique 
with one or more coherent states, as possibility of hav- 
ing more than N photons in the system potentially alters 
the postselection reliability. A starting point for such a 
study would be the detection model in Ref. [2ll |. 

There are various important considerations for the 
practical application of the postselected 2A^-port scheme 
we have presented. The first is that it can be built for 
TV < 4 with current technology. Down conversion sources 
capable of producing 4 single photons are becoming com- 
mon and production rates are increasing ,22i|; number 
counting photon detectors are also becoming available 
j2Qj. For small N , it is practical to replace the number- 
resolving detectors with beamsplitters and ordinary non- 
resolving photon counters, with slightly reduced overall 
probability of a postselection count. Taking iV = 3, for 
example, the number-counting detectors can be replaced 
with a non-resolving detector on output and a 50/50 
beamsplitter followed by two such detectors on output 
1. The postselection occurs on a triple coincidence, and 
from Eq. (|10|l it follows that the probability of such a 
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detection still varies with (j> as cos{N(j)). 

We have already mentioned the commercial availability 
of 8-port and higher fibre splitters. This, however, leads 
to an important consideration. For > 4, the general 
2Af-port device is not necessaril y d escribed by the canon- 
ical multiport matrix, Eq. (jH [ll|23|. For < 3, the 
conservation of energy requirement defines the 2iV-port 
matrix to within external phases on the input and out- 
puts, i.e. the canonical matrix entirely represents the 
physics for iV < 3 However, when > 4, there ex- 
ist free internal phases independent of the conservation 
of energy. These free phases lead to substantive changes 
in the 2iV-port matrix. As an example, the matrix for 
= 4, with free internal phase 9, can be written: 
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which reduces to the canonical case, Eq. ||SJ), when = 
tt/2. The existence of these free internal phases for A^ > 4 
provides an extra degree of freedom for making a wide 
range of nonclassical states using multiport beamsplit- 
ters with single photon and coherent state inputs. An 
open experimental question is how close the transfer ma- 
tricies of fused fibre multiports come to the canonical 
matrix, and whether production can be tuned to obtain 
the canonical case. 

Finally, there exists significant freedom to explore 
other interesting nonclassical states that may be pro- 
duced by techniques similar to the protocol we have pre- 
sented. We have already mentioned varying the internal 
phases of the 2A^-port device, as well as the use of more 
than one coherent state input, but one can also modify 
the output network, the postselection protocol, or the in- 
put states |23|. An instance of the latter is using n > 1 
Fock states. For example if the input state of an 8-port is 
set to = |2, 2, 1, 1), and single photons are detected 
in the bg and 62 modes, then in the modes 61 and 63 the 
state is exactly the NOON state, \iP)b = |4,0) + |0,4), 
without the need for postselection or zero photon detec- 
tions. This success probability of this configuration is 
3/64, the same as the more complicated circuit in Ref. 

In conclusion, we have theoretically demonstrated a 
technique for producing maximally entangled photon 
number states using 2A^-port fiber splitters, photon- 



counting postselection or conditional measurement, and 
Fock-state inputs. We have also shown that it is possible 
to replace one of the inputs with a coherent state and 
achieve the same maximally entangled state. This tech- 
nique allows a highly nonlinear state transformation to 
be realized experimentally with a simple optical circuit. 

Note added. Recently we became aware of an indepen- 
dent proposal for implementing phase measurements and 
nonlinear gates using multiport devices and postselection 
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APPENDIX 

We prove the result in Eq. (jSJ, with complex numbers 
(3, 7 replacing SJ, b\, since the latter commute. 
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Proof. Define the set of numbers Xk — (3 
These numbers are the roots of 
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which is in fact the characteristic equation for the matrix 



M = 
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(3 7 
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Thus we see that the left hand side of Eq. ljA.l|) is in 
fact the product of the eigenvalues for M, in other words 
it is the determinant of M. The determinant of M is 
simply the left hand side of the characteristic equation, 
Eq. l|X2|) , with A = 0. Thus Eq. (|Jt|i follows. 
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